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Abstract. We present a generahzation of continuous position measurements that 
accounts for a spatiahy inhomogeneous measurement strength. This describes many 
real measurement scenarios, in which the rate at which information is extracted about 
position has itself a spatial profile, and includes measurements that detect if a particle 
has crossed from one region into another. We show that such measurements can 
be described, in their averaged behavior, as stochastically fluctuating potentials of 
vanishing time average. Reasonable constraints restrict the form of the measurement 
to have degenerate outcomes, which tend to drive the system to spatial superposition 
states. We present the results of quantum-trajectory simulations for measurements 
with a step-function profile (a "which-way" measurement) and a Gaussian profile. We 
find that the particle can coherently reflect from the measurement region in both cases, 
despite the stochastic nature of the measurement back-action. In addition, we explore 
the connection to the quantum Zeno effect, where we find that the reflection probability 
tends to unity as the measurement strength increases. Finally, we discuss two physical 
realizations of a spatially varying position measurement using atoms. 
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1. Introduction 

Position measurements strike at the heart of what distinguishes quantum mechanics 
from classical mechanics. Within classical mechanics, it is possible to localize particles 
with arbitrarily fine precision for all times. In quantum mechanics, the back-action from 
the measurement disturbs the particle and implies that there are fundamental limits to 
measurement precision — limits that are being approached in current experiments [ , ] . 

The usual formalism of projective measurements fails when applied to position 
measurements: if a particle is reduced to an eigenstate of position it has an infinite 
momentum uncertainty and hence infinite energy. Position measurements are best 
described by considering the "weak" limit where the observer continuously extracts 
information about the particle's location at a given rate; in a small time interval only a 
small amount of information is extracted. 

Continuous measurements for quantum systems are usually phrased in the language 
of stochastic master equations, where the system's density matrix is updated based 
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on random measurement outcomes [ ]. The stochastic master equation can also be 
derived from a physical model for the interaction of the system with the measurement 
apparatus [ ]. The first model specifically for a continuous position measurement 
accounted for the ensemble- averaged evolution of the density matrix [ ] . This has been 
developed further into the modern description in which stochastic master equations 
(quantum trajectories) are used to continuously condition the quantum state based on 
the output of a classical detector that is correlated with the system's position [ ]. The 
observer can equivalently monitor the status of a large bath with which the system 
interacts. The scenario in which the resonance fiuorescence of an atom is monitored to 
gather position information is one example of this [7, 8]. A similar approach is used 
to treat measurements of the position of an atom interacting with the field mode of an 
optical cavity [ ] . 

Continuous measurements have an important role in quantum control and in 
the transition to the classical limit of quantum mechanics. Continuous position 
measurements in particular can be employed in feedback control loops for cooling 
quantum systems [ ]. This has been applied to atoms in cavities [12-15], trapped 
ions [16, 17], and nanomechanical resonators [ ]. In addition to seeking to control 
quantum systems, continuous position measurements also provide one path through the 
quantum-to-classical transition. In particular, continuously monitored quantum systems 
are able to exhibit chaotic behavior [19, 20], in contrast to closed quantum systems. 

Considering the importance of position measurements, it is necessary to develop 
the theory to account for realistic constraints applicable to any experimental realization 
of a position measurement. In particular, in any real position measurement the particle 
can only be detected within a limited region. If the particle leaves this region, the 
observer gains no further knowledge of the particle's position. This can be modeled 
by a space-dependent coupling of the particle to a bath, such as the radiation field. 
When the strength of the coupling to the bath, and thus the measurement strength^ is 
itself a function of the position, it is, in fact, this function of the position that becomes 
the measured observable. The coupling correlates the bath with that function of the 
particle's position, and so monitoring the bath allows the observer to gather information 
about that function. 

The average dynamics of the measured particle (that is, the motion of the particle 
averaged over all the possible measurement results) can be reproduced by a fiuctuating 
potential with the same position dependence as the measurement. Thus, as far as 
the average motion is concerned, the measurement acts like a stochastic force. This 
fluctuating potential provides an alternate, intuitive physical picture for understanding 
the measurement back-action. We expand on this unitary "unraveling" of the average 
motion in Sec. 2.2. 

It has been shown previously by a number of authors that a measurement that 
determines whether a particle is on one side of a dividing line or the other can exclude 
the wave- function from the region in which the particle is initially absent [^1-29]. This 
causes the particle to reflect from the dividing line, and is due to the quantum Zeno 
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effect. In our case this corresponds to a measurement with the foUowing two properties: 

1) the measurement strength takes one of just two values, and makes a jump from one 
value to the other at the dividing line (the measurement strength is a step function); 

2) the measurement strength (equivalently, the size of the step) is much greater than 
the kinetic energy of the incoming particle, suitably scaled. In this case the bath exerts 
a large back-action on the particle, providing the momentum transfer required for the 
reflection. While this Zeno effect has been discussed previously, we examine it here 
to connect it with the stochastic master equation formalism. We do this by deriving 
an analytic expression for reflection probability in Sec. 4. We also perform numerical 
simulations of a measurement with a step-function proflle, for comparison to a localized 
measurement with a Gaussian proflle. 

A position measurement can be physically realized by monitoring the resonance 
fluorescence of an atom interacting with a resonant light fleld. The proflle of the 
light becomes the measurement function, and the measurement strength is proportional 
to the overall intensity. As a consequence of the aforementioned effects, the atom 
will coherently reflect from the resonant light fleld for a sufliciently large intensity, 
despite the absence of a mean dipole force, and despite the stochastic nature of 
resonant atom-light interactions that tend to heat the atom. This phenomenon, besides 
being somewhat counterintuitive, has implications in situations where atoms encounter 
localized, resonant, optical-pumping flelds — as occurs, for example, in implementations 
of one-way barriers for atoms [30-3?]. We present further details on this realization in 
Sec. 5.1, along with another example of an atom interacting with an off-resonant cavity. 

We begin in the next section by deriving the stochastic master equation that 
describes a continuous position measurement with a spatially varying measurement 
strength, and elucidate some of its key properties. In Sec. 3 we perform simulations 
of a particle incident on measurements with two kinds of spatial proflles, showing the 
behavior of the particle on individual realizations of the measurement, which may involve 
either reflection or transmission, as well as the ensemble- averaged behavior. In Sec. 4 
we discuss the quantum Zeno effect, in Sec. 5.1 we present two physical realizations of 
a spatially varying measurement, and in Sees. 6 and 7 we flnish with some concluding 
remarks. 

2. Equations of Motion 

2.1. Derivation 

In this section we will derive the equation of motion describing a spatially varying 
position measurement. Our measurements are usually made by monitoring a large bath 
that interacts irreversibly with the system. The position measurement arises from a 
spatially-dependent potential that couples the particle to the bath. The bath then 
becomes correlated with a real function of the particle's position //(x) and we can 
distinguish different positions by how strongly they interact with the bath. Although 
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p — - — = p — t[H^ p] dt — KT-L[p\x)p{x)]pdt^ (3) 



p is real, we will keep our notation general because later we will add a complex local- 
oscillator amplitude to p. We will start our derivation from the positive-operator- valued 
measure for a measurement with two outcomes in each infinitesimal time-step, dt (note 
that our treatment and superoperator notation parallels that of Ref. [34]). We will refer 
to the two outcomes respectively as "no detection," and "detection" . (We choose this 
nomenclature due to the formal similarity to photodetection, though these outcomes 
are not necessarily tied to a photodetector.) The measurement operators that describe 
the two outcomes are 

= 1 — jH dt — Kp\x)p{x) dt (1) 
= V2i^dtp{x), (2) 

where H = p'^ /{2m) + V{x) is the free Hamiltonian for the system, k is the measurement 
strength, corresponds to "no detection", and to a "detection". The detection 
outcome occurs with probability Tr[J~^|f^ip] = 2k{p\x)p{x)) dt. The evolution that 
occurs between detections is due to f^o^ 

Ti[^l^^p] ~^~h' 
where the Ti super-operator is defined as 

n[c\p = cp + pc^ - {c + c^)p. (4) 

If detection occurs the density matrix changes according to 

^ ^ip^^l ^ p{x)pp\x) 
^ Ti[^\^^p\ Ti[p^{x)p{x)p\' ^' 

We can represent all of this evolution in a single stochastic master equation (SME): 

dp = —-[H^ p] dt — i<iH[p\x) p{x)]p dt + Q[p{x)]p dN ^ (6) 

where 

Here, dN is a Poisson process [35, 36], where dN = 1 with probability 2k{p\x)p{x)) dt 
and is zero otherwise. If we mix in a "local oscillator" with this signal, then we pass 
over to a master equation similar in form to the usual position-measurement master 
equation [ ] . This corresponds to the transformation 

a . . 

\/2k 



H^H ^(o^> - ap), (86) 

where a = l^le^"^ is the complex amplitude of the local oscillator, which leaves 
the unconditioned evolution (that is, the evolution averaged over all the possible 
measurement results) unchanged. For large we can pass over to the white noise limit. 
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When the detection rate becomes fast on the time scale of the particle's motion, the 
central limit theorem allows us to replace dN by a mean drift and Gaussian fluctuations, 

where double angle brackets denote the ensemble average for the random variable, 
and dW is Ito white noise, with {{dW)) = and dW'^ = dt [35, 36]. We apply the 
transformations in Eqs. (8a), (86), and the white noise limit of Eq. (9), to the jump 
master equation (6). In the large \a\ limit, the resulting white noise master equation is 

dp = -^[H, p] dt + 2K.V[p{x)]pdt + V2KrhL[e~'^p{x)]pdW, (10) 

where the V superoperator is deflned as 

V[c\p = cpc^ — ^ (c^cp + pc^c) . (11) 

For = 0, Eq. (10) is simply the standard stochastic master equation that describes the 
measurement of the observable O = p{x) [ ]. Thus the case p{x) = x recovers the usual 
master equation for a measurement of position. The continuous stream of measurement 
results, often referred to as the measurement record^ is given by 

dW 

dr = {p{x))dt + (12) 
\/8k 



where dr is the result in the time interval dt. 

The unconditioned dynamics are given by the master equation 

i 

dtp = -j^[H,p] - i^[p{x), [/x(x),p]], (13) 
where we have taken an ensemble average over all possible noise realizations. 



2.2. Stochastic Potential 

An alternative interpretation is suggested if we choose a different phase for the local 
oscillator. If we instead choose = 7r/2, the master equation becomes 

dp = -^[H dt + V2Khp{x)dW, p] +2KV[p{x)]pdt, (14) 

= -^[Hdt + V2ilhp{x) o dW, p], (15) 

where the flrst equation is in Ito form, and in the second equation o dW denotes 
Stratonovich white noise [3o, 3()]. This phase choice corresponds to measuring the other 
bath quadrature [37], yielding information about the measurement backaction rather 
than position. The effect of the bath here is formally equivalent to a stochastic potential. 
The spatial dependence p{x) of the coupling to the bath determines the proflle of the 
stochastic potential. The term in V[p{x)\ is the Ito correction to the Stratonovich 
fluctuating potential. When averaged over all noise realizations, this master equation 
recovers the unconditioned evolution given by Eq. (13). 
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This view is particularly apparent if we look at the input-output formulation of a 
homodyne measurement [ ]. The measurement arises from a linear coupling between 
the observable and the bath, 

i^SB dt = + dBl), (16) 

where dBi^ is a white noise bath operator due to the sum over bath modes, all evolving 
at different frequencies [39]. This interaction term correlates the state of the bath with 
The key observation is that since the interaction Hamiltonian is proportional 
to the noise in the bath drives the system through The bath therefore 

exerts a stochastic force on the particle whenever fi{x) is inhomogeneous. Note that 
any measurement of that is mediated by a bath must have an interaction of this 
form [ ]. This in turn implies that there must be a stochastic force of the same form. 
As a concrete example, a two-level atom interacting with a resonant laser field feels the 
stochastic dipole force due to spontaneous emission [il]. 

2.3. Effects of Measurement 

We can gain insights into the dynamics induced by a spatially varying position 
measurement by examining the lowest few terms in the Taylor expansion of about 
the current mean position {x). If is constant then the post-measurement state is 
completely unaffected, and the measurement does nothing. 

The linear term in the Taylor expansion acts as a standard position measurement. 
The measurement can then be approximated by 

V2i^lJi{x) ^ V2i^iJi' {{x)){x - (x)), (17) 

which defines the effective local measurement strength, where ii'{x) = dji/dx. The 
constant (x) term again has no effect. This linear approximation is valid if the particle 
is well-localized on the scale of the measurement function. This relation provides the 
link between a general, inhomogeneous measurement of position and the usual, linear 
position measurement. 

The quadratic term in the Taylor expansion is important at the maxima and minima 
of the measurement function where the linear term vanishes. The effective measurement 
is 

V2i^fi{x) ^ V^ilji" {{x)){x - (x))^ (18) 

which drives the state to spatial superpositions since the measurement result cannot 
distinguish {x) + x from (x) — x [ ]. Once the components of the superposition move 
away from the extremum, the local measurement for each component becomes linear. 
Each component will then evolve under their respective effective measurements. 

2.4' Form of the Measurement Function 

We now derive the conditions that the measurement function /i(x) must satisfy for it to 
correspond to a position measurement that acts only in a limited region of space. First 
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Figure 1. Possible measurement functions realizing a position measurement within a 
bounded region. The local measurement strength is proportional to the derivative of 
(a) Measurement function with different asymptotic values, (b) Measurement 
function with a discontinuous return to zero, (c) Continuous measurement function 
with same asymptotic values. 



we expand the measurement terms in the master equation (10) in the position basis, 
with p{x^x') = {x\p\x^): 

dp{x^ x') = — — x^) dt 

+ V2i^ [p{x) + /x(x') - 2(/x)] x') dW. (19) 

The leading term leads to decay of coherence between positions x and x' where p differ. 
If we enforce p{x) = /x(x'), then the measurement cannot distinguish between x and x\ 
This implies that /x(x) should be constant outside the measurement region. 

Inside the measurement region, the measurement function /x(x) should be linear 
to act as a standard position measurement. However, if /x(x) does not return to the 
same value on either side of its linear section (as in Fig. 1(a)), then the measurement 
would collapse superpositions that are on either side of, and completely outside the 
linear measurement region. This is not what one normally means when restricting a 
measurement to a particular region. For example, a photodetector of some finite size 
placed well within the arms of an optical interferometer should not collapse the fringes by 
giving which- way information, but this is precisely the case for the function in Fig. 1(a). 
Thus, for a measurement to truly act only within some bounded region, the function 
must take on the same value everywhere outside that region. We may take this value 
to be zero, since shifting fi{x) by a constant has no effect on the dynamics. 

Generally speaking, /x(x) will be the profile of some continuous (e.g., laser) field. 
Thus, a discontinuous return of the measurement function to the external value, as in 
Fig. 1(b), is at best an idealization of a more physical, continuous measurement function. 
Along the same lines, the measurement should also be finite in extent. Measurement 
functions of the type in Fig. 1(a) must therefore also be idealizations of measurement 
functions that are eventually zero, though possibly only returning to zero far away from 
the location of the measured particle. 

The strictly physical form of the measurement function must therefore be 
continuous and return to the same constant value on both sides of the measurement 
region, as shown in Fig. 1(c). In particular, this means that there will always be 
indistinguishable positions Xi, X2, where /x(xi) = p{x2) within the measurement 
region, corresponding to the same measurement outcome. At any given instant, the 
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measurement cannot distinguish between the positions Xi, X2, even if the measurement 
is not symmetric on either side of the extremum, for example if /x'(xi) 7^ /x'(x2). The 
measurement wiU thus create a superposition. However, if we also account for the 
Hamiltonian evolution then the backaction from an asymmetric measurement can give 
the components of the superposition differing momenta, and hence different positions at 
later times. Hence over an extended period of time the measurement could distinguish 
between the two positions. Likewise, we have argued that for a single measurement 
channel, a strictly physical measurement function must result in positions that are 
indistinguishable. However, multiple position-measurement channels — each with their 
own ambiguities — may be combined to remove all the ambiguities. These are examples 
of combining multiple measurements to gain information that cannot be gleaned from 
either measurement separately. A more familiar example is that information from 
position measurements at two different times may be combined to obtain information 
about momentum. 

2.5. Mean Momentum and Kinetic Energy 

We can derive equations of motion for the two lowest-order moments of the momentum 
probability density based on Eq. (10). If we assume the particle Hamiltonian is given 

by 

H=^+V(.). (20) 

we obtain 

d{p) = - {d^V) dt + y2^[(/xj9 + pfi) - 2(/x) {p)] dW (21) 
d{p^) = - {pd^V + d^Vp) dt + 2h\{{d^fif) dt 

+ y2^[(/x/ + //x) - 2{fi){p^)] dW. (22) 

In an ensemble average d{{p)) = —{dxV) dt^ so the average momentum is not changed by 
the measurement. For a free-particle ensemble, the average momentum is conserved — 
during a single trajectory, the measurement can exert a force, but the force vanishes 
on average. The measurement back-action enters as diffusion in the ensemble-averaged 
moment evolution, 

d{{p^)) = 2h^K{id,^if)dt, (23) 

with momentum diffusion coefficient = 2h^n{{dxlJiY) - This reflects that it is the 
gradient in ii[x) that gives relative position information. 

3. Simulations 

We now present some quantum-trajectory simulations that reveal some of the novel 
dynamics that occur under inhomogeneous position measurements. In the simulations, 
the particle starts in a minimum-uncertainty, Gaussian state with mean position 
(x) = Xo, mean momentum {p) = po and root-mean-square width ax- It is incident 
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on a measurement function from outside the measurement region, from the left- 
hand side. Since the atom starts in a pure state we initiaUy have maximal knowledge 
about the atom. Our results focus on the case when po/cTp ^ 1, which corresponds to 
a well-defined momentum. In the limit of small the particle simply passes through 
the measurement region with mild heating, so we also restrict ourselves to the regime 
of strong measurements, tv > pl/{2mh). 

Our simulations use a split-operator Fourier method to propagate the stochastic 
Schrodinger equation associated with Eq. (10). We have set h = m = 1 throughout our 
presentation here. We visualize the evolution of the particle under the measurement 
using the Wigner function for the particle [ to] . 

3.1. Step function 

We will first consider a step measurement function 



which is simple and analytically tractable. It also affords a comparison with a familiar 
example from textbook quantum mechanics: the potential step. This corresponds to 
the limit of an arbitrarily large gradient for a continuous measurement. This limit is 
unphysical, strictly speaking, but it is useful as an idealization of a strong measurement 
over a small region. 

Fig. 2 shows the temporal evolution of a wave packet incident on a static potential 
step V{x) = Vo^(x), with no measurement, as a context for understanding the 
measurement-driven evolution. The evolution can be described by the Schrodinger 
equation, or by setting = in the measurement master equation (10), and including 
the potential V{x). The height of the potential barrier in this case is smaller than the 
kinetic energy of the wave packet, so the particle can transmit over the barrier. The wave 
packet splits into two pieces, refiected and transmitted. The refiected wave packet has 
exactly the opposite momentum of the inbound wave packet. However, the transmitted 
wave packet is slower, since the particle lost energy climbing the potential step. The 
striated red and blue regions denote interference. The fringes between the refiected 
and transmitted wave packets show that they represent a coherent superposition. Note 
that the fringes between the two coherent components are oriented along the direction 
between them. This corresponds to the fact the fringes would yield an interference 
pattern in the appropriate marginal probability distributions [ ]. Missing fringes 
indicate a classical mixture of the two states. 

Fig. 3 shows the evolution of a wave packet incident on a measurement step 
//(x) = 9{x)^ but for a free particle with no explicit potential. The step measurement 
distinguishes whether the particle is on the left- or the right-hand side of the origin, but 
it does not resolve different positions on a given side. This means that superpositions 
confined to one side are unaffected by the measurement, while superpositions on different 
sides of the origin will collapse. This measurement is analogous to a "which- way" 
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Figure 2. Wigner functions for a particle incident on a static potential step with 
height Vq = 0.5, and without measurement. The wave packet has initial momentum 
Po = 1^ width (Ja; = 5, and initial position xq = —Scix. Time is measured in units 
of dx/pO' Ked values are positive, blue values are negative, and black is zero. The 
corresponding animation is included in the supplementary data. 
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Figure 3. Wigner functions for a particle incident on a step-function measurement 
profile jii{x) = 0{x) with = 5; the profile is the same as the potential in Fig. 2. The 
trajectories are generated with a white-noise-measurement unraveling, as given by the 
Eq. (10). The initial conditions are the same as in Fig. 2. (a) Single trajectory with 
reflecting measurement outcome, (b) Single trajectory with transmitting measurement 
outcome, (c) Ensemble average over measurement realizations. Red values are positive, 
blue values are negative, and black is zero. The corresponding animations are included 
in the supplementary data. 



Reflection of a Particle from a Quantum Measurement 



12 



measurement in that it has only two outcomes: left or right. The two possible 
measurement outcomes correspond to reflection or transmission of the wave packet. 

Fig. 3(a) shows a single trajectory where the wave packet reflects from the 
measurement. Initially, we know the particle is to the left of and propagating towards 
the step in the measurement function. As the wave packet evolves, we never measure 
it to be on the other side, and this measurement outcome effectively forces the wave 
packet to reflect. The momentum transfer for the reflection comes from the back-action 
of the measurement. As noted in Sec. 2.2, a spatially dependent measurement implies 
a spatially varying back-action on the system. 

Fig. 3(b) shows an example where the particle transmits through the measurement 
step. The wave packet starts to reflect, but still has a small evanescent tail that 
penetrates the x > region. The interference between the reflected and initial wave 
packets, before transmission occurs, can be seen in Fig. 4. Random measurement 
outcomes consistent with the particle being in the x > region are rare at flrst since the 
evanescent tail represents a small probability. Once they occur, however, these outcomes 
can increase the amplitude of that part of the wave function. This can lead to a runaway 
process where future measurement results are more likely to indicate the particle is in the 
X > region, and thereby collapse the particle entirely into this region. This process 
occurs continuously on a time scale /^~^. Since the measurement time scale is faster 
than that of the motion in the strong-measurement limit, the particle is localized in the 
small, evanescent tail before it can propagate away. The particle will now have a broad, 
positive momentum distribution and will propagate away from the step, as shown after 
transmission occurs in Fig. 4. The transmission process is incoherent: the state of the 
particle is entangled with the state of the bath, which destroys any interference with 
the incident wave packet. Since transmission occurs at random times, the transmitted 
wave packets will have different phases and propagate away at different times. 

Fig. 3(c) shows an ensemble average of 10^ trajectories for particles incident on 
the step measurement. Note that the specular component of the reflection survives the 
ensemble average, while other features related to the measurement noise are washed out 
in the average. The transmitted parts show no coherence with one another or with the 
reflected parts. The incoherent parts of both reflected and transmitted wave packets 
show heating from the measurement back-action. 

Surprisingly, part of the reflection is coherent. It is possible to see interference 
between the reflected and incident wave packets, even though the reflection is caused by 
the back-action of a measurement, which conventionally leads to diffusion or heating. 
The interference is marked by the alternating positive and negative regions centered at 
p = in the Wigner function. The interference fringes have the same phase for each 
member of the ensemble — that is, they are the same for each possible, random outcome 
of the measurement — producing the coherent reflection. 

The measurement-induced reflection is analogous to interaction-free measure- 
ments [ ]. In those experiments, light detects the presence of an absorber in the arms 
of a polarization interferometer without ever interacting with the absorber. An optical 
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Figure 4. Position and momentum distributions for a particle before and after it 
transmits through the step measurement. The initial conditions are the same as in 
Fig. 2. The time elapsed between "before" and "after" is At = 0.4/^~^. Note that this 
is not the same trajectory as in Fig. 3(b), but it displays a similar behavior. Inset: log 
plot of the same momentum distribution. 



example of an interaction-free measurement closely related to the reflection phenomenon 
here is that of a high-flnesse Fabry-Perot cavity [45, 46]. An empty cavity transmits 
light almost perfectly, but when an absorber is placed in the cavity to "detect" the hght, 
a large fraction of the light instead reflects from the cavity, although the reflected part 
has not interacted with the absorber. 



3.2. Gaussian measurement 



We now consider the case of a Gaussian measurement function. This is an even function 
of position, so we expect this measurement to tend to drive the particle towards spatial 
superposition states. In addition, if we consider the resonance-fluorescence scenario, 
a laser probe with a Gaussian proflle would realize this Gaussian measurement. The 
measurement function is 

.2 



exp 



X 



(25) 
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where sets the width of the measurement. Note that we take to be dimensionless, 
which is important to keep in mind when comparing to standard measurements = 
X, which makes sense only when x is taken to be dimensionless. 

Fig. 5 shows a wave packet incident on a Gaussian potential barrier without any 
measurement. The wave packet has enough potential energy to cross the barrier (without 
tunneling). As with the step potential, the wave packet coherently splits into reflected 
and transmitted components. Note that in this case, both parts of the wave packet have 
the same asymptotic momentum after scattering from the potential barrier. 

Fig. 6(a) shows a single trajectory where the majority of the wave packet reflects 
from the Gaussian measurement, without an explicit potential. Since the Gaussian is a 
smooth function, the wave packet readily enters regions where the measurement strength 
is small. The Gaussian function lacks a sharp edge, and it is possible to get a series 
of results that weakly localize the atom in the soft edges of the function. A small part 
transmits through on this particular trajectory, as indicated by the striations showing 
the coherence between the reflected and transmitted parts of the wave packet. 

Fig. 6(b) shows an alternative case where the wave packet is localized by the 
Gaussian measurement. The sides of the Gaussian act as a linear measurement for 
wave packets conflned to one side. Each side can distinguish different positions within 
that side and act to localize the wave packet inside the measurement region. The wave 
packet then diflFuses around from the momentum kicks of the measurement. This has 
occurred by t = 4 in Fig. 6(b). If the wave packet passes over the peak of the Gaussian 
then it will be driven to a superposition state, since the measurement is locally even 
in X at the maximum. If the wave packet is moving quickly then the superposition 
will tend to be biased in the direction of motion. The creation of a superposition state 
is probabilistic, and depends both on localizing the wave packet, and the wave packet 
diffusing over the peak of the measurement function. After the superposition forms, the 
two components diffuse away from the peak of /x(x). The symmetry of the measurement 
ensures that the superposition is not disturbed once it is created. A superposition state 
has been created by t = 6 in Fig. 6(b). 

Fig. 6(c) shows the ensemble average over 10^ trajectories for a Gaussian 
measurement. Once again, the coherent reflection survives the ensemble average. The 
interference fringes between the reflected and incident wave packets are easier to see, 
since the incoherent features from weak localization and partial transmission have 
averaged out. Both the reflected and transmitted parts of the wave packet have an 
incoherent, heated part; this comes from realizations where the particle is localized and 
split into a superposition state. In general, the phase of each superposition is random 
and the superpositions are created at random times, so the coherent features of each 
superposition are washed out in the ensemble average. The ensemble- averaged results 
do not condition on any measurement scheme, and correspond to the unconditioned 
density operator — the solution to the unconditioned master equation (13). Thus we 
will see coherent reflection from any process that realizes master equation (13) in the 
ensemble average. 
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Figure 5. Wigner functions for particle incident on a static Gaussian potential, with 
potential width = 1 and height Vq = 0.25, without measurement. The wave packet 
has initial momentum po = 1^ width (Ja^ = 5, initial position xq = —3{(Jx + ct/^)- Time 
is measured in units of {cix -\-crfi)/po- Red values are positive, blue values are negative, 
and black is zero. The corresponding animation is included in the supplementary data. 
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Figure 6. Wigner functions for particle incident on a Gaussian measurement with 
measurement width = 1 and strength = 5; the profile is the same as the potential 
in Fig. 5. The initial conditions are the same as in Fig. 5. The trajectories are 
generated with a white-noise-measurement unraveling, as given by the Eq. (10). (a) 
Single trajectory with reflection from measurement, (b) Single trajectory producing 
a spatial superposition state, (c) Ensemble average of measurement. Red values are 
positive, blue values are negative, and black is zero. The corresponding animations are 
included in the supplementary data. 
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3.3. Heating 

We can estimate the momentum transferred by the measurement by considering the 
Fourier transform of the measurement function. This is related to the transverse 
momentum distribution of the bath particles that interact with the particle undergoing 
the position measurement. In the resonance-fluorescence example, the bath particles 
are photons that form the mode proflle of the intensity distribution which becomes 
the measurement function We obtain the rate of momentum disturbance by 

multiplying the Fourier transform of fi{x) by the measurement rate which is the 
rate at which we gather position information — and hence disturb momentum. The 
approximate flnal momentum distribution will be the convolution of the particle's 
momentum distribution with the Fourier transform of the measurement function. 
For example, the Fourier transform of the step function 9{x) is 



where P.V. denotes the Cauchy principal value. The sharp edge of the step implies 
an arbitrarily flne resolution of how close the particle is to the edge, which implies 
a divergent momentum disturbance. The convolution with 1/k leads to long tails in 
momentum space and a divergent momentum uncertainty. 

The convolution of momentum distributions as a result of "which-way" 
measurements has been noted before in the context of projective measurements [47] 
and transverse momentum transfer. The convolution result is exact for a projective 
measurement, but must incorporate dynamical evolution for continuous measurement. 

The Fourier transform of a Gaussian is also Gaussian, so the momentum transfer 
in this case is flnite. The width of the momentum distribution scales as the inverse of 
the width of the measurement resolution. A narrow measurement function constitutes 
a locally stronger measurement of position and hence a larger disturbance. 

There is an alternative, measurement-theoretic explanation for the heating. A 
strong measurement ensures that we detect if the particle crosses between regions of 
differing Consider the time when the particle begins to cross into the next region, 

as in the t = 2 and t = 4 panels of Fig. 3(b), or the t = 2 panel of Fig. 6(b). If the 
measurement record does not yet indicate that the particle has crossed through, then 
the particle can only have penetrated a small distance into the next region. Then, when 
the measurement indicates that the particle has passed into the next region, the state 
collapses down to the small part that penetrated the next region. Since this part is 
tightly conflned in position, there now is a large uncertainty in momentum. The wave 
function is effectively multiplied by /x(x), and so the wave function will acquire the 
same character in momentum space. This is apparent, for example, in the t = 4 frame 




(26) 



of Fig. 6(b). 



Reflection of a Particle from a Quantum Measurement 



17 



1 




03 

O 











50 



100 



150 



200 



Figure 7. Measurement-induced reflection probabilities for a wave packet with ax = 
10 incident on fi{x) = 6{x). Numerical simulations (averaged over 256 trajectories) are 
compared with the analytical result from Eq. (32). 



4. Quantum Zeno Effect 

Despite the heating, a particle will almost certainly reflect from a strong measurement. 
As Fig. 7 shows, the probability for reflection tends to unity as the measurement strength 
increases. This is a manifestation of the quantum Zeno effect [21-24, 48, 44, 49]. 

A simple argument in the case of a step- function measurement = 9{x)^ as 

in Sec. 3.1, demonstrating this in the limit of a sequence of projective measurements 
separated by infinitesimal time intervals (/^ — > oc) is as follows ['^^]. We assume that 
at time t the particle is completely confined to the region x < 0, as enforced by the 
projective measurement in the immediate past. The operator for evolution between 
measurements is 



The Hamiltonian is responsible for any evolution that could transfer the wave packet 
into the region x > 0. The term proportional to dt is the probability amplitude for 
the particle to cross into the region x > during the infinitesimal time step, while the 
amplitude for staying in the x < region is unity. The probability of detecting the 
particle in the region x > is thus 0{dt^) = for the next projective measurement, 
resetting the particle to the x < region and thus confining it there. A similar argument 
has been made before from the perspective of a "gap" in quantum mechanics [ ] . 

A more quantitative argument can be used in the case of finite measurement 
strength. This is related to quantum reflection of atoms dropped at grazing incidence 
onto ridged surfaces [51, 52]. This "quantum reflection" has been related to the Zeno 
effect since the ridges act heuristically as absorbers for the atoms, and thus increase 
the reflection probability from the surface [ ] . This treatment misses the measurement 
back-action and heating that we consider here. Our calculation here also makes a 
rigorous connection of this absorber model to continuous-measurement theory. 




(27) 
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Here we will model the evolution by a jump-process realization for = ^(x), 

2 

d\^) = - dt-K {e\x) - (e'^ix))) dt 

+ (J}MT^-W)<iN, (28) 



where dN is a Poisson process with mean 

{{dN)) =2K{e\x))dt. (29) 

This stochastic Schrodinger equation gives the same unconditioned evolution as the 
white-noise unraveling in Eq. (10). We will model the incoming wave packet as a 
unit-amplitude plane wave, which gives the correct reflection coeflicient as long as the 
momentum distribution is narrow compared to the mean incident momentum. 

We can use a linear trajectory and solve for the between-jump evolution (the 
nonlinear term represents the probability of the particle to be detected in the region 
X > 0, which is precisely what we are computing here). The linear, between-jump 
evolution is governed by 

^2 



ihdf 



ihnOix) 



(30) 



We assume the particle is initially located on the left with mean energy pl/2m^ and 
enforce continuity of the wave to match the incoming wave with a reflected wave of the 
same energy and a decaying solution for x > 0, which gives the total probability for the 
particle to be in the "transmitted" region. This is used along with the mean detection 
rate to flnd the total detection probability: 



Pdet = 2/^ j dx9{x)\ij{x)\^ (31) 

— oo 

= 272^/^ .32) 

where ^ := 2mhn/p^ is the measurement strength scaled in units of the incident kinetic 



energy, and X •= y ^1 + — 1. For large the detection probability, and thus the 
transmission probability, vanishes as . 

In this linear-trajectory picture, the measurement can be viewed as an imaginary 
potential that absorbs particles. The absorbed particles are detected inside the x > 
region, and their wave functions are modifled according to 

1^) ^ ^^"^1^^ (33) 

The remainder of the particles coherently reflect from the measurement with a 
probability given by the norm of the state. 

In Fig. 7 we have counted only the coherent part of the reflection by computing 

J dp |(/>(— +]^)p5 where ^(p) is the momentum-space wave function. This omits the 

— 3(7t, 
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contribution from the diffusely reflected parts. As the initial momentum increases, the 
particle penetrates farther into the x > region, and hence the transmission probability 
increases. 

This calculation works because the step measurement is a two-outcome 
measurement — as soon as the particle is detected it collapses into the measurement 
region. This is analogous to the calculations that can be made for the probability of 
spontaneous emission from a two-level atom [ ] . This approach can in principle be used 
to calculate the coherent reflection probabilities given that the particle remains outside 
the measurement region for a general measurement function. 

Although the reflection probability can approach unity, the ensemble-averaged 
momentum is conserved, as shown in Eq. (21). The conservation of mean momentum 
ties together the observation of increasing reflection and the heating of the transmitted 
wave as increases. As tv increases the probability for reflection increases as the wave 
packet is excluded from the measurement region. Concomitantly, any increase in 
increases the heating from the measurement since the evanescent tail in the absorbing 
region is more tightly conflned. This oflFsets the decrease in mean momentum from the 
reflection. 

The Zeno eflFect has also been noted when combining a continuous position 
measurement of = x with a double- well potential [ ]. For moderate measurement 
strengths, the position measurement localizes the particle in one well and inhibits the 
tunneling between the two wells. This breaks down for large measurement strengths 
since the heating from the measurement boosts the energy of particle above the barrier 
in the middle. In our case the Zeno eflFect arises because of the inhomogeneity of the 
measurement function, rather than the inhomogeneity of the potential. 



5. Physical Realizations 

5.1. Stochastic Dipole Force 

The above effects can be realized by coupling a two-level atom to a high-intensity 
light fleld with position-dependent Rabi frequency. This can be viewed in terms of 
the stochastic dipole force [ ]. We write the position-dependent Rabi frequency as 
Q{x) = \Q\g{x)^ where \Q\ = \deg • E|//i, |E| is the maximum electric fleld amplitude, 
and g{x) is a complex, dimensionless mode proflle with \g{x)\ < 1. The atom obeys 



dtp = -i 



+ AcrV + - [n*{x)a + n{x)a''] , p 



+ rV[a]p, (34) 



_2mh 2 

where F is the excited state decay rate. In the high intensity hmit, the dressed states 
|±) are approximately equal superpositions of the ground state |g) and the excited state 
|e), and vice versa. Each spontaneous-emission event projects the atom into the ground 
state, and thus into an equal superposition of the dressed states. The atom therefore 
sees both shifts ±h\fl{x)\/2 of the respective dressed states |±). The interpretation is 
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as follows: after each spontaneous emission event, the atom sees a force 

F ^ ^5^. (35) 

where the sign is chosen randomly, but with equal probability for the two possibilities. 
Assuming the atom is moving slowly, the momentum change associated with a single 
spontaneous-emission event is 

Ap^T?^« (36) 

where ^ is the time until the next spontaneous emission event, which is a random variable 
> 0) of mean 2/r and exponential probability density 

/(0 = ^exp(^-^^). (37) 

To take into account the randomness of the sign, we can write 

Ap^'^f (38) 

where G M has a two-sided exponential probability density 

/±(0 = 5exp(^-^iri). (39) 
Then the mean-square kick is 

«^,).) = Mr) = ™^. (40) 

where 

The diffusion rate is the mean-square step divided by the average step time At = 2/r, 
so 

^„^<^f^^M^q^i«'(.)p. (42) 

This heuristic result agrees with the high-intensity limit of more rigorous but 
cumbersome calculations [^^^ 41^ F ]. 

Since the mean time r/2 between changes in the force is much smaller than 
typical time scales for the motion of cold atoms, we can think of the atomic motion 
in a temporally coarse-grained sense. Then the stochastic evolution of the atomic 
motion over a time interval of interest involves many transitions of the potential. The 
momentum changes have finite variance, so we can invoke the central limit theorem, 
and treat the force fiuctuations as (Gaussian) white noise. An equivalent, Stratonovich 
stochastic potential that gives rise to the same diffusion rate is thus 

V{x,t)=^-^V\9Wo^, (43) 

where we have dropped the mean dipole potential, which vanishes on resonance A = 0. 
We can see that this form of the potential is correct by comparing Eqs. (14) and (15), 
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so that the Ito form of the SME has the same form as Eq. (14), with = |r2p/2r and 
= ^y g^{x)g{x)^ where the sign of is defined such that \fi\x)\'^ = \g\x)\'^. We 

can then read oflF the momentum-diflFusion coefficient using Eq. (23) to see that it agrees 

exactly with Eq. (42). 

It is interesting to note that in this high-intensity regime, the stochastic potential 

represents an interplay between the ffuctuations in the atomic dipole and the field 

intensity. That is, this effect is due to ffuctuations in the atom-ffeld interaction 

Hamiltonian itself. As such, it does not saturate for high intensities, as does spontaneous 

emission. 

A measurement interpretation of the stochastic dipole force is as follows. The 
dipole force for an atom in a particular dressed state arises due to coherent scattering of 
photons between different wave vectors of the driving ffeld. In principle, this coherent 
redistribution of light can be measured to obtain information about the position of the 
atom in the ffeld, as for example is the case in recoil-induced resonances [5^^]. However, 
on resonance, the atom ffips randomly between the dressed states, so that the mean 
redistribution of photons (and hence the mean dipole force) vanishes. The measurement 
information is instead encoded in the variance of the redistributed light. 

Note that a low-intensity jump decomposition of resonance ffuorescence has been 
considered in the context of realizing imaginary potentials [ ] , though without exploring 
the consequences for the atomic dynamics. Unfortunately, in the low-intensity limit the 
spontaneous emission recoil obscures the measurement effects, and this model can break 
down when pushed to the high intensities required for a large measurement strength. 

5.2. Cavity Quantum Electrodynamics 

A spatially varying position measurement can also be realized in the interaction of a 
two-level atom with the ffeld of an optical cavity, which is in turn driven by a classical 
ffeld. Due to a large detuning, and a large cavity decay rate we can eliminate the 
atom's internal dynamics and the cavity dynamics, leaving the only the center-of-mass 
motion of the atom. In this case the cavity mode function becomes the measurement 
function [9, 13, 14]. The simulations from Sec. 3.2 correspond to Gaussian mode of a 
ring cavity. 

We assume that the cavity is being monitored by homodyne detection, so the system 
obeys the following master equation [l4] 

dp = -^[H, p] dt + rV[a]pdt + jV[a]pdt + ^/jn[a]pdW, (44) 

where F is the free-space spontaneous emission rate, 7 is the cavity decay rate, and a is 
the cavity annihilation operator. The total Hamiltonian is given by 

H = ^ + /2iJca^a+ncjA^V+%(x)(aV+aat) + n£;(ae^^^'+ate-^^^'), (45) 

where uJc is the cavity mode resonance, is the atomic transition frequency, ^'(x) is 
the cavity mode proffle, and E = yJ'^Pjhjj^ is the amplitude for the classical driving 
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field, where P is tlie power in tlie driving field. If we transform to a rotating frame, 
we can remove the free evolution due to the cavity. Since the classical driving field is 
resonant with the cavity, we can then write 

H = — + hAa^a + %(x)(aV + aa^) + hE{a + a^), (46) 
2m 

where A = — c^c- In the limit A ^ ^'(x), 7, F, p^/(2m/i), we can adiabatically 
eliminate the excited state. This amounts to replacing a and with (a) and {cr^). (For 
a more rigorous approach see Ref. [ ].) In this limit the effective Hamiltonian is given 
by 

//eff = — + h^-^a^a + hE(a + a^), (47) 
2m A 

and the spontaneous emission terms can be dropped since the effective decay rate is 
0[A~^]. If we further assume that the cavity is strongly damped such that 7 ^ 
^'^/A, p^/(2m/z), we can also eliminate the cavity. The final result is the following 
master equation for the center-of-mass dynamics. 



where 



dp = -■^[//eff,p] dt + 2/^P[/x(x)]prft + y2^7{[/x(x)] dW, (48) 



Ii{x) 



2 4 

9l 



(49) 



where a = 2£'/7, and 5^0 = max[5'(x)]. The mean potential in the effective Hamiltonian 
can be cancelled by the Stark shift from an off-resonant classical field that does not 
resonate with a cavity. This interaction has the form 

^^Stark = (50) 

Cancelling the mean potential allows the measurement effects we have discussed to come 
to the fore. The master equation then has the same form as Eq. (10). 



5.3. Imaginary Potentials 

Measurements can also be viewed as imaginary (absorbing) potentials, as seen in Sec. 4. 
Imaginary potentials have been used to coherently diffract atomic beams [59, 60] and 
reshape wave functions to counteract the expansion of a free particle [ ]. A resonant 
standing wave of light acts like an imaginary potential for two-level atoms. When the 
atom absorbs a photon from the resonant standing wave, it will decay to another state, 
thus being effectively absorbed from the initial beam. The diffraction and reshaping 
both rely on post-selecting atoms that have not undergone spontaneous emission. The 
surviving atoms are thus most likely to be near the nodes of the standing wave. For 
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an atomic beam, the standing wave acts as a diffraction grating and enables Bragg 
diffraction for atomic beams. This so-caUed "anomalous diffraction" has been further 
explained in terms of degeneracies of the imaginary Hamiltonian [b2, ]. The shaping 
can be viewed as conditioning on null detection of a sin^(/cx) measurement of position, 
while the diffraction is a consequence of scattering off the periodic array of nodes. 

6. Analogies 

The coherent reflection from a measurement is analogous to the coherent back-scatter 
of light from a disordered atomic medium [ ]. The analogy is clearer from the 
fluctuating-potential picture, since the paraxial wave equation is the same as the 
Schrodinger equation with time replaced by the propagation direction. Under the same 
correspondence, the fluctuating potential in time becomes a disordered potential in 
space. In the ensemble average the fluctuating potential also yields a coherent reflection, 
since it obeys the same unconditioned equations as the measurement. 

An alternative analogy is the reflection of a photon from a conducting surface. The 
large complex permittivity — whose imaginary part represents absorption — leads to a 
rapid extinction of the wave inside the conductor if the conductivity is large. This is 
analogous to a large measurement strength for detecting the electromagnetic wave inside 
the region. This leads to the large reflection probability for light from the surface of a 
good conductor, even though it is a good absorber for waves inside the medium. 

The above examples capture the coherence of the reflection, but miss the 
transmission aspects of the measurement since the inbound particle is absorbed by the 
interaction. Unlike photons, atoms are not destroyed upon detection, so the form of the 
transmitted wave packet is necessary for a complete picture of measurement-induced 
reflection. 

7. Outlook 

We have outlined a theory of continuous position measurements that describes a spatially 
varying measurement strength. The effects shown here require the measurement gradient 
to dominate over all other dynamics. It is only in the limit of a large measurement 
gradient on the scale of the atom, large intensities and small velocities that coherent 
reflection can take place. This work can be extended to include other couplings and 
position measurements such as atoms coupled to micro-toroidal resonators [ ] , and the 
effect of measurements when the particle is confined to a potential. The case of both an 
inhomogeneous measurement and potential is of interest, for example, for studying the 
quantum-classical transition for classically chaotic systems. Future work will extend 
the idealized theory presented here to account for the imaging setup, diffraction effects 
and realistic photo-detection. 
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